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Abstract. - The DC current- voltage characteristics induced by a driving electric field with 
frequency Q of a one dimensional electron channel with a tunnel barrier is calculated. Electron- 
electron interaction of finite-range is taken into account. For intermediate interaction strengths, 
the non- linear differential conductance shows cusp-like minima at bias voltages mfiQ/e (m 
integer) that are a consequence of the finite non-zero range of the interaction but are independent 
of the shape of the driving electric field. However, the frequency-scaling of the photo-induced 
current shows a cross-over between and Q,~ 2 , and depends on the spatial shape of the 
driving field and the range of the interaction. 



The influence of a mono-chromatic driving field on the transport through tunnel barriers 
has become an extensively studied subject. A paradigmatic result has been obtained more 
than three decades ago by Tien and Gordon for a superconducting tunnel junction Similar 
approaches have been used to study photo-assisted effects in the tunneling through quantum 
point contacts ^ and driven quantum wells || . Basically, the frequency of the driving field has 
been shown to produce sidebands due to the non-linearity of the current-voltage characteristics. 
In all of these works, interaction between the electrons has been neglected. Recently, the 
influence of electron-electron interaction on non-linear tunneling transport has been studied in 
the regions of charging, interaction Q, ||, || and strong correlations ||, [h]] . In mean field 
approximation, renormalisation of the driving field has been studied |llf| . Also experimentally, 
in view of the realization of the quantum dot current turnstile device |jL2[ , photo-induced 
tunneling through quantum dots has been often addressed |ll| |lj, [l^, [ll| . 

In this paper, we present physically striking results from evaluating the general theory 
developed earlier, applied to driven non-linear transport through a tunnel barrier in a one 
dimensional (ID) system of interacting electrons. We report a novel frequency-locking effect, 

Typeset using EURO-TeX 



2 



EUROPHYSICS LETTERS 



which is signature of the systems coherent, strongly correlated electron states. It is character- 
istic of the finite, non-zero range of the interaction but does not depend on the exact shape of 
the driving electric field. 

By starting from the assumption that the current is driven by a local electric field of 
arbitrary shape, with a DC-component E (x) and an additional mono-chromatic component 
with frequency Cl, E(x, t) = E a (x) + E\(x) cos Clt, we find an extremely rich behaviour of the 
photo-induced current-voltage characteristic I(Vq; V\,Cl) as a function of various parameters, 
like the range and strength of the interaction and the spatial shape of the applied field (Vq, 
Vi voltages corresponding to Eq, E\, respectively), (i) We first re-derive for finite range 
interaction the non-trivial general expression for I(Vq; V%, CI) in terms of the (non- linear) DC 
current-voltage characteristic Iq(Vq) = I(Vq; Vi = 0, Cl — 0) || [l?], [is) 

I(V ;V 1 ,Cl)= J2 JZ(\z\)I (y a + n™y (1) 
n=—oo ^ ' 

It generalises the Tien-Gordon result to strongly correlated tunneling objects and arbitrary 
shapes of the driving fields, (ii) We show the existence of pronounced, cusp-like minima in 
the differential conductance d//dVb at integer ratios of eVo/HCl if the dimensionless interaction 
parameter g > 2/3, independent of how exactly the driving voltage drops. These structures are 
due to sidebands that are induced by the non-linearity of the DC current- voltage characteristics 
[ [l7[ jlj}] , and are independent of the shape of the driving electric field. Finally, (iii) we 
predict a cross-over between frequency scaling cx f2 _1 and cx Cl~ 2 for delta function-like and 
spatially constant driving fields, respectively. In the cross-over region the frequency scaling is 
non-universal and reflects the ranges of both, the electron-electron interaction and the range 
of electric field. 

It is well known that interaction-induced renormalisation of the externally applied field 
influences the current in the frequency and time domains [^[ [1^, |ll| . Our results indicate that 
for ID correlated electrons, as long as the driving field is mono-chromatic, the Tien-Gordon 
formula remains valid, but with a general argument of the Bessel functions. The latter is given 
by the modulus of the complex quantity 



/oo 
AxE 1 {x)r{x,Cl), 



and contains the spatial shape of the AC-componcnt of the driving field, and the interaction 
range via the AC-conductivity a(x, Cl) without tunnel barrier, r(x,Ct) = a(x,Cl)/a(0,Cl). 

In order to study the effect of finite range interactions and a space-dependent electric field on 
non- linear AC-transport, we use the ID Luttinger liquid model JO], ^fj. The Hamiltonian 
for a spin-less Luttinger liquid with an impurity and subject to a time-dependent electric field 
is H = Hq + H t + H ac , where Ho is the Hamiltonian describing charge density excitations 
with a dispersion uj{k) — v-p\k\[l 4- V ee (k) /Httvf} 1 ^ 2 (vf Fermi velocity). It reflects the Fourier 
transform of the interaction potential p3] . We assume a 3D screened Coulomb interaction of 
range a -1 projected onto a quantum wire of the diameter d. The dimensionless interaction 
parameter of the model is g = [1 + V cc (k = 0) /hnvF}^ 1 ^ 2 ■ For a -1 ^d, the interaction decays 
exponentially and one obtains the "Luttinger limit", V cc (x) cx ae""' 1 ' p^ |. 

The tunneling barrier of height Ut is assumed to be localised at x — 0, 

H t = U t cos [2Vtti?(x = 0)] (3) 



with the phase variable $(x) giving the charge density fluctuations in the limit of long wave- 
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lengths p(x) — k-p/ir + d x -d(x)/ y/ir. The coupling to the driving field is 



H ac = e 



dxp(x)V{x, t). 



(4) 



The voltage V(x,t) is related to the field by E(x,t) = —d x V(x,t). For the space-dependence 
of the electric field we assume Ei(x) = E\e~^ x ^ a , with the voltage V\ = — dxEi (x) = 
—2E\a. When the range of the field tends to zero, E\(x) reduces to a ^-function. The spatial 
dependence of the DC part of the electric field does not need to be specified, as only the voltage 
Vq = — dxE (x) is of importance in DC transport 

The current at the barrier is given by the expectation value I(x = 0,t) = (j(x = 0,i)), 
where the current operator is defined via the continuity equation, d x j{x, t) = —edtp{x, t). For 
a high barrier, the tunneling contribution to the current can be expressed in terms of forward 
and backward scattering rates which are proportional to the tunneling probability A 2 . The 
latter may be obtained in terms of the barrier height Ut by using the instanton approximation 
p2| . The result can be written in terms of the one-electron propagator S + \R Jul , 



I(x = 0, t) = eA 2 J dre- s ^ sin R{t) sin | £ 



dt'V eS (t') 



with 
S(r) 



ii? M = ^ / < ^-Ke{a- 1 (x = 0,u;)} 



{ 1 — cos ut) coth — 



lsmuT 



(5) 



(G) 



where (3 = l/k^T, w max the usual frequency cutoff that corresponds roughly to the Fermi 
energy |23|| , and the AC conductivity of the system without impurity is |9j 



a(x, ui) 



hir 2 



°° cos kx 

'^(^-(u + io+y 



Furthermore, the effective driving voltage is related to the electric field by (|] 



v eS (t) 



dx 



dt'E{x,t')r(x 1 t - t') =V Q + — \z\ cos (fit - ip z ) , 



(7) 



(8) 



where \z\ and Lp z are, respectively, modulus and argument of z (cf. Eq. (|^)). The modulus \z\ 
will be the argument of the Bessel functions in the final result, (p z will represent a phase shift 
of the harmonics. With the above assumptions about the shapes of the driving field and the 
interaction potential one obtains 



eVi 



-.A 



n k(n) 



y/l + a 2 k 2 {ny \v F a 
where k(Q) is the inverse of the dispersion relation and 



-, aa 



A 2 (u, v, w) = 



1 



1 



, (u + wv) 2 
(uw + v) 2 



(9) 



(10) 



In the following, we concentrate on the results for the DC component of the current which 
does not depend on x and is directly given by the current at the barrier, for which we only 
need to know Izl, 



= eA : 



: / dre~ s,(r) sin R(t) sin 
Jo 



eVnr 



Jo 



2|z|sm— J 
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Fig. 1. - Currents To, Idc and differential conductance d/dc/dVb at zero temperature as a function of 
the ratio eVo/hQ for g = 0.9 (left), g — 0.5 (right) for values fl — v^a, a = 0, and eV\/hv^a — I 
(£ = 5 dotted, I = 6 dashed, 1 = 1 dash-dotted lines). Currents in units of hvFa/eR t ; differential 
conductance in units of ; tunneling resistance Rt = 21100^^/ ne 2 A 2 . 



(11) 



n— — oo 



The important point here is that the driven DC current is completely given in terms of 
To(Vb), the nonlinear DC current-voltage characteristic of the tunnel barrier, 

I (Vq) = eA 2 j dre- s(r) sin R{r) sin (^rj . (12) 

Eqs. (|Tl|) , ( p^| ) generalise results which have been obtained earlier |Q but without interaction 
between the tunneling objects, and also for the Luttinger model with a zero-range interaction, 
together with a 5-function like driving electric field || . 

For Vq much smaller than some cutoff-voltage V c which is related to the inverse of the 
interaction range, I oc V* /9 -\ This recovers the result obtained earlier for ^-function 
interaction and zero-range bias electric field (?[. When V » V c , the current becomes linear 
[|8|. For intermediate values of Vo, lo exhibits a cross-over between the asymptotic regimes 
with a point of inflection near V c . For zero-range interaction, Iq oc V^ 9 1 for any Vq. For the 
understanding of the driven current-voltage characteristic to be discussed below, this general 
shape of Iq(Vo) for a finite-range interaction will turn out to be very important. 

Figure [l] shows the currents Iq, Idc and the differential conductance didc/d^o as functions 
of eVo/hil for g = 0.9 and g — 0.5 for zero-range of the driving electric field. For g — 0.9 
one observes sharp minima in the differential conductance at integer multiples of the driving 
frequency in certain regions of the driving voltage V\ . These can be understood by the following 
argument. When the strength of the interaction is not too large, the region where didc/dVo 
is much smaller than 1 is small compared with ft£t such that for eVb ~ didc/dVo oc 
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Fig. 2. - Scaling exponent v of the argument \z\ of the Bessel functions as a function of Q, for ranges 
of the driving field (curves from right to left) aa = 10~ 3 , 10~ 2 , 10"\ 1, 10, 10 2 , 10 3 , and g = 0.5. 



(2/g - 1) |eVo - W\ 2 ^- 2 . Then, Eq. (|T]) yields near eV = mKl 

^ « 1 - + const • Jl{\z\) \eV - mh^ 2 . (13) 

For g > 2/3, this yields for integer m the cusp-like structures observed in Fig. [I]. For g < 2/3, 
no cusps occur anymore. In addition, the current Jd c is depleted so strongly and over such a 
large region of the bias voltages that the regime of almost vanishing d/dc/dVb becomes larger 
than Ml and in general no minima near integer multiples of the frequency exist. As can be seen 
in the figure, the depths of the cusps depend on the driving voltage V\ (oc \z\) which can also 
be understood from of Eq. ( |l3| ) which shows that the values of the differential conductances 
at the voltages eV — mtel are approximately 1 — J^(|z|). 

It is therefore instructive to look into the behaviour of \z\ as a function of the frequency. 
Figure || shows the scaling exponent v determined from 

dloglzl , „ 

dlog$2 

We observe a non- universal cross-over between \z\ oc Q^ 1 , the case discussed by Tien and 
Gordon which corresponds to a driving field of zero-range (a — > 0), and \z\ oc Q,~ 2 which is 
obtained for a homogeneous external field (a — > oo) H. Although the behaviour of z depends 
strongly on the parameters of the model in the cross-over regime, this does not influence 
qualitatively the occurrence of the cusps. Their existence depends crucially on the finite range 
of the interaction, and the condition g > 2/3. However, by varying \z\, the depths of the 
minima are changed due to the variation of J^(|.z|). 

Given the above result for the driven DC-current, the general behaviour of the differential 
conductance as a function of eVo/hfl can be straightforwardly obtained. Of special interest is 
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the occurrence of cusps at eVo/Ml = m (m integer) which appear to be quite stable against 
changes in the model parameters. A similar result has been discussed earlier ]2j|, but for a 
small potential barrier between fractional quantum Hall edge states which implies zero-range 
interaction. In the general case discussed here, the finite range of the interaction is crucial for 
obtaining the cusps, due to the absence of a linear contribution towards the current for small 
voltage which is characteristic of tunneling in ID dominated by interaction. The cusps could 
be used to frequency-lock the DC part of the driving voltage. 

In conclusion, we have demonstrated that the result which has been obtained by Tien and 
Gordon for tunneling of non-interacting quantum objects in ID driven by a mono-chromatic 
field localised at the tunnel barrier remains valid even in the presence of interactions of 
arbitrary range and shape, and for an arbitrary shape of the mono-chromatic driving field. The 
central point is that the frequency driven current is completely given by a linear superposition 
of the current-voltage characteristics at integer multiples of the driving frequency, weighted 
by Bessel functions. 

The argument of the latter contains the amplitude of the driving voltage only linearly but the 
dependence of the argument on the frequency and the range of the driving field is determined 
by its spatial shape. However, one can easily identify regions where the dependence on the 
frequency becomes very simple. For a driving field which is localised near the tunnel barrier, the 
integral in Eq. (|^) can be evaluated approximately by noting that r(x, Q) varies only slowly with 
x and can be taken out of the integral. Then, \z\ — eVi/Ml which corresponds to the result of 
Tien and Gordon |l]]. In the other limit of an almost homogeneous electric field, Ei = Vi/a, one 
needs to calculate the spatial average of r(x, f2) Q. This gives a(k = 0, Q)/a(x — 0, f2) » fi -1 , 
since a(x = 0,O) » const. This implies \z\ oc f2 -2 . Such a frequency dependence has been 
discussed earlier for non-interacting particles Here, we see that it is valid under quite 
general assumptions also for interacting particles. A possible method to detect this behaviour 
experimentally is to investigate the real part of the first harmonic of the current through the 
tunnel contact and to determine the current responsivity which is given by the ratio of the 
expansions of and the first harmonic to second and first order in \z\, respectively pp| . 

Finally, we wish to comment on the question of self-consistency Jll|] . The driving field has 
been assumed ad hoc in the present paper, and not determined self-consistently by taking 
into account the internal interaction-induced re-arrangement of charges and currents as a 
consequence of an external field. There are basically two effects that have to be considered. 

First, re-arrangement of the charges and currents will lead to a change of the spatial shape 
of the driving field. As discussed above, the present results allow to identify features that 
are independent of the spatial shape of the field and thus can be expected to survive even if 
self-consistent corrections are taken into account. Second, apart from the displacement current 
contribution, which is present even for A — > and renormalises \ z\, additional contributions to 
the local field with other frequencies may be generated, especially in the present nonlinear case. 
Even the DC current could be changed due to mixing of harmonics present in the local field. If 
self-consistency maps contributions with different frequencies into the local field, possibly the 
above result have to be generalised. Strictly speaking, the Tien-Gordon form will no longer 
be valid. However, as has been discussed earlier [[To), the amplitudes of the nth harmonics 
generated by the tunnel barrier when starting from a mono-chromatic external voltage decay 
extremely rapidly with n. Thus, when using them as a starting point for local field corrections, 
practically only the fundamental frequency would contribute. This would eventually lead only 
to changes in the spatial behaviour of E\ and not destroy the existence of the cusps. In 
summary, we claim that local field corrections which need to be treated in any case beyond 
the mean field approximation are in the present context unimportant and would not lead to 
essential changes of the our findings. This concerns especially the existence of the cusps in the 
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differential conductance at very low temperatures that are characteristic of the non-linearity 
induced by the finite-range interaction. But further, more detailed studies are necessary, in 
order to clarify the quantitative aspects of this issue [p8| . 
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